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I. INTRODUCTION 

In this work we derive the equations of motion for spin- 
ning test bodies in the context of gravitational theories 
with nonminimal coupling. Some of these theories have 
recently been investigated in [IH3- I n particular, we gen- 
eralize and extend the findings in Q to the case in which 
the nonminimal coupling depends generally on the curva- 
ture of spacetime. Our results apply to a very large class 
of theories, i.e. the coupling is allowed to be a general 
function of the set of 9 parity-even curvature invariants. 
The covariant multipolar expansion method used to de- 
rive the equations of motion in the present work goes back 
to @ - it has also been utilized to derive the equations of 
motion of test bodies in Einstein's theory [6] . Our find- 
ings therefore allow for a direct comparison of the new 
'force' terms due to the nonminimal coupling procedure 
in a covariant fashion. 

Without going into historical detail, we only men- 
tion that in the context of General Relativity multipo- 
lar methods of the kind employed in this work - and 
variations of it - were used in the works [6 15]. Similar 
methods have also been successfully applied in alterna- 
tive gravity theories, see [l6l - fl8| and more recently in 
(l9l [20l | . A short timeline of works can be found in [19[ . 

Note that the model under consideration does not be- 
long to the very general class of gravitational models ana- 
lyzed in [19J due to its nonminimal coupling prescription. 

The structure of the work is as follows: In section HI1 we 
introduce two nonminimal coupling scenarios; in particu- 
lar, we work out their respective conservation laws. With 
the conservation laws at hand, we derive the equations 



of motion for test bodies - up to the pole-dipole order 
- by means of a covariant multipolar method in section 
IIII1 This is followed by our conclusion and outlook in 
HVi Appendix [X] contains a brief overview of our conven- 
tions and notation, the dimensions of all quantities in the 
paper and a directory of used symbols can be found in 
tables U and ITT] Appendix [B] summarizes the expansion 
formulas used in our derivations. 



II. THE MODELS UNDER CONSIDERATION 

A. Nonminimal f(R) gravity 

In [l[ an extended version of a so-called f{R) grav- 
ity theory was considered. Gravity theories in which the 
usual Einstein-Hilbert Lagrangian is replaced by an ar- 
bitrary function of the curvature scalar have attracted a 
lot of attention during the last few years see, e.g., the re- 
views pll - [23| and references therein. The f(R) scenario 
was generalized even further in [l[ by the introduction of 
a nonminimal coupling term on the Lagrangian level. In 
particular, the following Lagrangian was put forward: 



fi(R) + [l + \f 2 (R)]L mat . 



(1) 
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Here /1 and f 2 are arbitrary functions of the curvature 
scalar R, and L mat is the matter Lagrangian. The non- 
minimal coupling of matter and gravity is controlled by 
the constant A. The general field equations - in terms of 
the functions /1 and f 2 and their derivatives - are given 
in their explicit form is irrelevant for the subsequent 
analysis though. 

In contrast to standard f(R) gravity theories, the last 
term in {T]) leads to a modification of the equations of 
motion. As was already shown in [24| (see also eq. (5) of 
fH) the usual conservation law - as, for example, found 



2 



in General Relativity - is replaced by 



V 4 T„ = 

1 1 + A/ a 



(gijL ma t — Tij) V l i?. 



(2) 



Here f 2 (R) '■= rf/2 (-R) /di? denotes a shortcut for deriva- 
tives of the unspecified function fa (R) of the curva- 
ture scalar and the energy-momentum tensor of mat- 
ter is defined in the standard way by ^J—gTij := 
-25(^L mat )/5g i3 . 



B. General nonminimal gravity 



The above model can be generalized to 

(3) 

where both the gravitational field Lagrangian L grav = 
L sia . v (gij,Rijk l ) and the function F = F(g ij ,R ijk l ) can 
depend arbitrarily on the spacetime metric and the Rie- 
mannian curvature tensor. For example, in Q both are 
assumed to be the functions of the Gauss-Bonnet scalar 



G = R A 



ARijR 10 



RijkiR 1 ^ 



(4) 



This case belongs to the general class of models when 
the Lagrangian L grav — L grav (ii, 13) and the function 
F = F(ii,i2 1 iz) both depend on the quadratic scalar 
invariants constructed from the components of the cur- 
vature tensor, 



ii = R\ 



1 2 



h = RijkiR v 



(5) 



As is well known |25l - l27j , a curved spacetime manifold of 
4 dimensions is characterized by the 14 algebraically in- 
dependent invariants constructed from the components 
of the Riemann tensor. There are two types of invari- 
ants: some of them are built of only the metric gij and 
the curvature Rij k l , while others involve also the Levi- 
Civita totally antisymmetric tensor. The invariants of 
the first type are parity-even quantities (i.e., they do 
not change under space and time reflections), whereas 
the second type of invariants are parity-odd objects that 
change their sign under coordinate transformations which 
do not preserve orientation. There are 9 parity-even in- 
variants and 5 parity-odd ones [26|, [2?J . Here we will 
confine our attention to the general nonminimal coupling 
theories in which F = F{i\, . . . , ig) is an arbitrary func- 
tion of the set of the parity-even invariants that includes, 
in addition to the quadratic scalars (JSJ) , the following cu- 
bic, quartic and quintic contractions: 

«4 = Rij kl Rki mn Rrnn 3 , (6) 

«5 = R l jR : 'kR k i, «6 = R l jR J kR k iR l i, (7) 
^ 7 = RVDij, i 8 /),,/)". i 9 (8) 

Here we have denoted := R ik ijR kl . The set ©-© 
is equivalent to the one reported in [2^, [27j, when the 



Riemann tensor is decomposed in terms of the Weyl and 
the traceless Ricci tensor. 

Generalized gravity theories with Lagrangians which 
are functions of the minimal independent set of curvature 
invariants have recently attracted some attention in the 
cosmological context, see [28|, for example. 

The field equations for the model © are derived from 
the variation of the total action with respect to the space- 
time metric. Denoting ^J—gEij := 2S(y/~ gL glav )/Sg t3 , 
we find explicitly 



E 13 = FT 13 
9 



+2 [L ma tF A P% + V„V fc (L mat F A n 



ikjn 
A 



(9) 



Here Fa — dF/diA, A = 1, . . . , 9. The second line de- 
scribes the modification of the gravitational field equa- 
tions due to the nonminimal coupling. Here, for the cur- 
vature quadratic invariants (A — 1,2,3), 



if = -R lk R 3 k -R lkl3 R ku 



p'j 

■*3 



-2RR" 3 , 

1 r>jklrn 
IK hlmft 



and 



ikjn 

ikjn 
T J — 



r> Ty{ in jk kn ji\ 

2R(g 5 - g g ), 

R in g jk _ R kn g ji _ R ij g nk 



R k3 \ 



(10) 

(11) 
(12) 

(13) 
(14) 



ir tk3n = 4R ik3n . 


(15) 


For the homogeneous cubic and quartic 


invariants 


«4, i 5 , i 6 , given by ©-(0), we find 




P? = -ZR lkln R ln ^R 3 kpq , 


(16) 


3 (2) "\ (2) 

Pi 3 = -~RhR ik -±R iMj Rki, 


(17) 


(3) (3) 

P^ = -2R 3 k R tk -2R mj R kll 


(18) 


(2).. ., . (3).. 

where we denoted R %3 :— R tK R 3 k and R 13 := 


K Ukirt J , 


TT lk3n = QR lk pq R 3npq , 


(19) 



r^'" = l(R in g 3k 



(2), .. 

R kn g 31 



(2 



■u- ( 3 )- ( 3 ), 

T Lk3n = 2(R m g 3k - R kn g 31 



R l3 g nk 



(3)-- , 
R V g nk 



(2), . . 

R k3 g m ), 



(20) 



R k °g m ) 



(21) 



Furthermore, for the mixed Riemann-Ricci quartic and 
quintic invariants © we derive 



Pi 3 = -D lk R 3 k - D 3k R l k - R lkl3 D ku 



( 2 ), 



Pi 3 ' = —Y ik R 3 k — Y 3k R i k — R ikl3 Y k i — D ij , 



(22) 
(23) 



P^ = -V lk R 3 \~V 3k R l k -R tkl3 V kl -^R m ^D k i 



^{R kl D\D 3 l + R l k D 3k 



R 3 kD lk ), 



(24) 
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(2) 

where F iJ := R lkh ' D kl and D 13 := D tk D^ k , and 



ikjn 
7T 7 


= R in R? k - 


- P fc "P^ 










+ D m g jk - 




D ij g nk ^ 


D kj g n \ 


(25) 


ikjn 
^8 


= D in RP k - 


- D fe "P Jl 


- D tj R nk 


+ D kj R ni 






+ y«yfc_ 


Y kn gji _ 


yij g nk + 


Y k3 g n \ 


(26) 


ikjn 
^9 


= U in R jk - 


- U kn R ji - 


- U zj R nk 


+ U kj R ni 






+ v*v*- 


ykn g ji _ 


yij g nk + 


ykj g ni 






1 (2). 


(2) (2) (2) 
k _ D kn g ji _ D ij g nk + D kj g 


ai \ 



(27) 



Here U ij := R^ k D'^ k and F y := R ikll U u . 

One can verify the symmetry properties P 
and 



ikjn 



ik]jr 



ik[jr 



jnik 



A ■ 



(28) 



for A = 1, . . . , 9. It is also straightforward to prove that 
the last term in (J9j) is symmetric in the indices 

The tensors ([ID ]) - ([15 ]) and (|I6 > ([2"? 1) satisfy certain dif- 
ferential identities. The latter arise from the fact that the 
action-type integrals I a = J d 4 Xy/—giA, A = 1,...,9, 
are invariant under general coordinate transformations. 
The corresponding Noether identities read 



\ - ( \iA9 ij 



A 



V n V fe 7T 



ikjn \ 



= o. 



(29) 



ViVnVfcTr^'" 



In view of the skew symmetry 

V[jV n ] VfeTr^'™. Then, using the definition of the cur- 
vature (IA1I). we rewrite the last term as 



ViV„Vfc7T* 



-Kikl V„7T A 



(30) 



As a result, the Noether identities ((29)) are recast into 



V7 . P % 3 



1 



1. 



^Rikl°^n^T n = ~ 2 V ' (31) 

Actually, one can verify these differential identities di- 
rectly by using the expressions (fTU|) - (TT51) and (|16I) - (|2"T)) 
for A = 1,...,9. 

The nonminimal model ([3]) is invariant under the 
diffeomorphism (general coordinate) transformations. 
The corresponding Noether identities then tell us that 
Vi-E 1 - 7 = identically, whereas on-shell (i.e., when the 
matter field equations are satisfied) 



V J FT i] + 2 



9 

£ 

A=l 



■^mat-PU-P^ 



r 771 ikjn 

^mati'A^A 



0. 



(32) 



We can simplify this considerably by making use (|3"0)) 
and the Noether identities (|3~Tj) . A direct check shows 
that flQ|) - (Pj) and (fl6 > p? l) satisfy (for A = 1, ... , 9) 

' ' ""■ (33) 



A 



/ > j_ klni 
-.n-kln K A 



0. 



Taking into account this crucial relation, the conservation 
law (l32t is recast into the final form 



V l Tjj — — (gijL ma t — Tij) V*P. 



(34) 



This result generalizes the conservation law © to the 
case in which F — F(i±, . . . , ig) depends arbitrarily on 
the complete set of 9 parity-even curvature invariants 
©-(IS]). Our derivation corrects the earlier studies @, @|- 



III. EQUATIONS OF MOTION 

In the following section, we derive the multipolar equa- 
tions of motion for test bodies from the conservation law 
(I34p by means of the covariant expansion technique from 
[5j. The multipolar moments extend the ones encoun- 
tered in [6J] to the general nonminimal coupling scenario. 

The equations of motion will be explicitly worked out 
up to the dipole order; i.e., they are applicable to general 
spinning test bodies in theories with nonminimal cou- 
pling. 



A. Covariant moments &: multipolar expansion 



To begin with, we rewrite (|3"4"|) as follows: 



V l 7V, 



ij - L/ 'mat 



Tn) V l A. 



(35) 



Here we have introduced a scalar function 
A[gij,Rijk l ) := log P. In the following discussion, 
we are going to denote derivatives of this function 
simply by Ai :— WiA, Aij := VjViA, etc. Raising the 
indices and rewriting the covariant derivative in (|35[) 
yields 



(36) 



In the last equation, we introduced the quantity 3 y := 
g lJ L ma ,t as a shortcut. Densities of different quantities 
are denoted by a tilde " " . 

We will now derive the equations of motion of a test 
body up to the dipole order by utilizing the covariant 
expansion method of Synge [5|. For this we need the 
following auxiliary formula for the absolute derivative 
of the integral of an arbitrary bitensor density B XlVl — 



D 

ds 



P* lw d£ T 



S(s) 



S( S ) 



•" ; v. v ,/>" v </ N :, . 



(37) 



E(s) 



Here v Vl := dx Vl /ds, s is the proper time, and the inte- 
gral is performed over a spatial hypersurface. Note that 



4 



in our notation the point to which the index of a bitensor momentum tensor T x ° Xl as follows: 
belongs can be directly read from the index itself; e.g., 
y n denote indices at the point y. Furthermore, we will 
now associate the point y with the world-line of the test 
body under consideration. For additional comments re- 
garding the explicit calculation of w a see the appendix 
of @. Now we start by integrating (|3"6")l : 



S( S ) 



a yi... a yn g yo xo lp : 



T x ° x AA Xl w X2 dZ x l38) 



V X1 (a y (/"•,, / ' ' „•'••,/>:, 



S( S ) 

S( S ) 

1 ,»»+ l r B ox 1Vw ((T m . ..^-gw.^)^. (39) 

S( S ) 

Equation (l39l) allows us to rewrite the integral (|38|) , and 



Here a denotes Synge's world-function and ff » its th f eb ^ to derive the e q uations of motion at arbitrary 



first covariant derivative, cf. also appendix [X] for a brief 
overview of our conventions. With the help of (|37[) we 
can rewrite the integral over the derivative of the energy- 



order. 

We now introduce integrated moments a la Dixon in 
i, i.e. 



n Vl---VnVO 



j.V2---Vn + iVayi — 



tV2---yn+iyoyi 



(-1)" / a yi ---a y »g y \ T x ^dX Xl , 

S( S ) 

(-1)" / * y2 ---<T y "+'g y ° X0 g y \ 1 T x ° x iw x *<E X2 , 

S(«) 

(-1)" J o y2 ---o y ^g y \ g y \ 1 % X0 ^w x *<E X2 . 

E( S ) 



(40) 
(41) 
(42) 



Then the equation ([38]) together with (|39|) , and the covariant expansions of the derivatives of the world- function and 
of the parallel propagator (see also appendix [B]) , yields 

_pVi—ynVo — j-(yi—y n )ya _ v (yipV2—yn)yo _ —R y ° y >yi>y tf»»n+u/V _)_ pVi---yny n +iy' v y"^j 

OO r 

_i_ \ , ( fyi---ynyn+i---y n +ky'yo -f-yi---y n y-n+i---y n +ky'yo\ a , 

+ / , £j ^? C J ^y'Vu+l'-Vn+k 

k=2 ' L 



n 3/2— »»)sn+i~ y n +kyo„,y' 



(43) 



In the dipole order, we keep only the multipole moments 
constructed up to the second order in the world-function 
cr. This truncates the infinite set of equations (|43p . with 
n = 0,1,..., oo, to the three lowest-order equations. 
Namely, for n = 2 and n = 1 we find 



v {yipy2)yo _ ^{yiV2)vo 
D 



(44) 



c/.s 



pj/i2/o _ ^j/oyi _ v yip 



+ Ay (£y i y°y 2 — t yiyoV2 ^ 



(45) 

The constraint equation (|44|) actually coincides with the 
one found in the general relativistic case. 
Finally, for n = 0, we obtain 



d« P ~ 2 y2!/l!/3 



t *3 



f^Vipysyz 

+ A yi (£ yoVl -F 1 ) 

+ A y y (j^yzyvy 1 — t V2VoVl ) 



(46) 
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B. Rewriting the equations of motion 



Equations (|44l) - (|46|) are all the information up to dipole 
order which we can extract from the integrated energy- 
momentum conservation law. The set of these three 
equations should be compared to the corresponding set 
of equations, which we derived in [i[ in the context of Pa- 
papetrou's [8j method. To make this even more explicit, 
we rewrite (pH |) - (|46p as follows: 

_a 1 r/a J.oi i A {tab +ab\ 

-j^P — 2^ bcdV s + A b[t; - 1 ) 



+ A bc (£ CQb - t cab ) , (47) 

—s ab = 2v [b p a] + 2A C ( t^ c - . (48) 

as V / 



Here we have introduced the spin of the test body under 
consideration as s ab :— 2p^ ab \ and switched back to the 
usual tensor notation, keeping in mind that all indices 
are now taken at the base-point y, which parametrizes 
the - still completely arbitrary - world-line. 

The equations of motion for the momentum (1471) and 
the spin (|48p should be compared to our previous find- 
ings (|3|,23) and ([H,29) in the context of a non-covariant 
multipole method. 

By utilizing the symmetries of the integrated t, p, and 
£ moments the dipole equations of motion can be further 
rewritten as 

E V " = lR\ cd v b S cd + C b ^bF + £ cab V c V b F,(49) 
as l 



—S ab = 2v [b V a] + 2£_ [ab]c V c F. 
ds 



(50) 



Here we have introduced the generalized momentum and 
spin tensors as 



■pa 
S ab 



Fp a +p ba V b F, 
Fs ab . 



(51) 
(52) 



Note that in (1491) and (|50|) the t moments have been com- 
pletely eliminated. 

A general interesting aspect of the present multipolar 
method is the fact, that the generalized momentum V a 
follows from the equations of motion. To recall, within 
the context of Papapetrou's method we had to specify 
this quantity by hand to achieve the final form of the 
first equation of motion. In the present approach it is 
retrieved from the equation of motion for the spin, i.e. 
([50]) yields: 



V a = Mv a + S ab v b - 2£ [ab]c v b \7 c F, 



(53) 



with M. := V a v a = Fm + p ab v b V a F ', where as usual, 
m := p a v a . 

It is worthwhile to mention the dual role played by the 
nonminimal function F: On the one hand, it "rescales" 
the ordinary momentum, spin and mass; and on the other 
hand, its gradients determine the force and torque that 



act on a particle in addition to the usual gravitational 
and Mathisson-Papapetrou forces. 

Finally, we note that - as expected - the motion of 
single-pole test bodies is also non-geodesic for the class 
of models under consideration. In this case, the geodesic 
equation, as encountered in General Relativity, is re- 
placed by 



£ (Fmv a ) = t ab V b F 
as 



(54) 



This can also be rewritten in an equivalent form 

mv a =^{5l-v a v b )V b A. (55) 

Here we recall that 3« := g ij L mat , and thus £ ab = g ab £„ 
with £ := J L ma ,tW X2 dY, X2 . As we see, a massive 
£(«) 

particle moves non-geodetically under the action of the 
"pressure" -type force (|55|) produced by the nonminimal 
coupling function F. 

IV. CONCLUSION 

In Q we employed Papapetrou's Q - non-covariant - 
approach to derive the equations of motion of the the- 
ory proposed in [JJ . The method utilized in the present 
work is more straightforward and has the benefit of being 
covariant. In 4 dimensions, there exist 14 algebraically 
independent curvature invariants [25j ]. The results ob- 
tained in (|49l) and (|50p generalize our previous findings 
to the case in which the nonminimal coupling depends 
arbitrarily on the Riemannian curvature of spacetime, 
with F = F(i\, . . . , ig) being any function of the com- 
plete set of 9 parity-even curvature invariants ([S])-©. 
The remarkably simple conservation law (|34p derived in 
this work corrects the earlier erroneous results [2, 3]. Fur- 
thermore, our final equations of motion explicitly make 
clear, that the previous non-covariant method relies on 
a very specific transport process; i.e., the choice of the 
parallel propagator g Vo Xo - It is satisfying to see that we 
formally obtain the same equations of motion with the 
more general method, which still allows for a recovery of 
the previous results in a special case. 

It is worthwhile to note that our results are compati- 
ble with the well-known general relativistic equations of 
motion for a spinning test body; i.e., in the minimal cou- 
pling case they reduce to the ones in 0, which in turn 
have also been derived by several authors by means of 
different multipolar approximation schemes 
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TABLE I: Dimensions of the quantities. 



Dimension (SI) 



Symbol 



TABLE II: Directory of symbols. 



Symbol 



Explanation 



Geometrical quantities 



m 

™2 



Matter quantities 



1 

kg/m 2 s 
kg m/s 
kgm 2 /s 
kgm n+1 /s 



Auxiliary quantities 



1 






-2 


nr 






-4 


rrT 






■6 


nT 


m~ 


■8 




10 


m~ 



kg/m 2 s 

_-n+2 



Operators 



s, dx a 
a 

Pabc , Pab, P 



rpab T J^ ao ^ a 

p a ,V a ,m,M,t 

„ab c<ab 



pab "z^ab 



ci... c n ab . . .c n ab ^Cj . . .c n a 



( 2 ), 



F, I A , A, A/ 2 

__abcd 
^1,2,3 

abed s~i ■ T)ti b 7~i 1 7 V, 

^4,5,7, tr, H,2,3, ^1,2,3) F> , R 

(3) 

abed pab ■ r rij -yij j-,ij 

(2) 

abed nab • \rii 

^9 , P 6,S, *6,8, V J , D 
T~>ab 

•Fj.,2,3 
#1,5,7 
F 6 ,8 
F 9 

h 

n V0 flVO ryVO 



Geometrical quantities 



gab 



x a , s 

Habc 

a 



I'D 



Metric 

Determinant of the metric 
Kronecker symbol 
Coordinates, proper time 
Connection 
Curvature 
World-function 
Parallel propagator 
Gauss-Bonnet scalar 



Matter quantities 



v 

m 

P a 

gal 
rp a I 



c\ . . .c n ab 



e 1 - 



Velocity 
Mass 

Generalized momentum 
Spin tensor 

Energy-momentum tensor 

Lagrangian 

Integrated moments 



Auxiliary quantities 



A 

h, h, F 

IA 

Fa 

Pa, K A bcd , U lj , V lj , 

.. (2).. (2).. (3).. 

D X1 , R' 3 , R 13 



Coupling constant 
Functions of the curvature 
Scalar curvature invariants 
Deriv. of F w.r.t. to invariants 
Shortcuts 



'yi-'-yn, P vo yi...y n , l Vo vt—y n Expansion coefficients 



Operators 

_D _«.„ 
ds 



(Partial, covariant) derivative 
Total derivative 
Coincidence limit 



Appendix A: Conventions & Symbols 



Our conventions for the Riemann curvature are as fol- 
lows: 

2T Cl - Cfc dl ... di . [H = 2V [a V fc] T Cl - Cfc dl ... d! 



' di...di 



" ^ ^ Rabdj T- k d\ . ..e...di ■ 



(Ai) 



The Ricci tensor is introduced by Rij = Rkij k , and the 
curvature scalar is R = g^Rij. The signature of the 
spacetime metric is assumed to be (+1, —1, —1, —1). 



In the following, we summarize some of the frequently 
used formulas in the context of the bitensor formalism 
(in particular for the world- function a(x, y)), see, e.g., p| 
129 . |30| for the corresponding derivations. Note that our 
curvature conventions differ from those in p, |30( • Indices 
attached to the world-function always denote covariant 
derivatives, at the given point, i.e. a y := V a cr, hence 
we do not make explicit use of the semicolon in case of 
the world-function. We start by stating, without proof, 
the following useful rule for a bitensor B with arbitrary 
indices at different points (here just denoted by dots): 



\B.. 



= \B.. 



\B.. 



(A2) 
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Here a coincidence limit of a bitcnsor B,,,(x, y) is a tensor 
[B.J = lim B...(x,y), (A3) 

x->3/ 

determined at y. Furthermore, we collect the following 
useful identities: 

a yoyix y 2 x 1 = cr yoyiy 2 x xi = a x xiy yiy 2 1 (A4) 

tf* 1 * 2 ^^ - 2a - ^^V yi a W2 , (A5) 
[a] - 0, [(7,] = [a,] - 0, (A6) 

[Cxia^] = [Cyjyj] = 9y 1 y 2 , (A7) 
[CaJijte] = Ps/iad = .9j/li/2 1 (A8) 
[Crr^jajg] = [c 2 ; l2:2 j )3 ] = [fTa^jfej/s] = [^J/i^to] = 0) 

(A9) 

[s x V] = [s x V;s 2 ] = = o, (A10) 



\o x ° 1 — -/? yo 

yy yi;x 2 x 3 j — r, 11 - yiy 2 y3- 



(All) 



Appendix B: Covariant expansions 

Here we briefly summarize the covariant expansions 
of the second derivative of the world-function, and the 
derivative of the parallel propagator: 



9 V *i (-<J«V 

OO s 



T yo 



2/1 



yi 

OO 



k=2 



yiy 2 ---yk+ 



ia y2 ---a Vk +\ (B2) 



n ya _ „v a y _r>y a y 3 

y xi;x 2 — y x t y x 2 \ y'y"yz u 



OO 1 v 



fe=2 



O yo — a y' ( - B Vo / rr y3 

il xi;y 2 — y xi I 2 J1, 2/ 2/2 2/3 u 



OO s 

"E ^7 y %'y 2 y3...// fc+2 ^ 3 ---^ fc+2 ).(B4) 
fe=2 ' ' 



The coefficients a, f3, 7 in these expansions are polynomi- 
als constructed from the Riemann curvature tensor and 
its covariant derivatives. The first coefficients read as 
follows: 



fyVO — _If?2/0 , 

" y\y 2 yz ~~ g Jl (2/23/3)2/1' 



y\y 2 yz 



1 



nyz)y\ 1 



yiy 2 ysy4 ~ ^ 



1 



(y 2 rL 2/32/4)2/1) 



«2/o — 7? yo 

I- 1 2/12/22/32/4 r, v (2/2-"' 2/3 3/4) 3/1' 



(S/S 71 |2/l|2/4)3/2- 



(B5) 

(B6) 

(B7) 
(B8) 

(B9) 



In addition, we also need the covariant expansion of a 
usual vector: 



00 f-1l fc 



fc=0 
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